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Abstract 

We study the spin-ordering and the magnon collective modes of the two- 
dimensional Wigner crystal state at strong magnetic fields. Our work is 
based on the Hartree-Fock approximation for the ground state and the time- 
dependent Hartree-Fock approximation for the collective modes. We find that 
the ground state is ferromagnetic, i.e., that all spins are aligned at T = even 
when the electronic g-factor is negligibly small. The magnon calculations show 
that the spin-stiffness is much smaller in the crystal state than in fluid states 

which occur at nearby Landau level filling factors. 
PACS numbers: 73.20.Mf, 75.30.Ds, 75.30.Et 



1 



Typeset using REVTgX 



I. INTRODUCTION 



At sufficiently low electron densities or in sufficiently strong magnetic fields, electrons 
will crystallize at low temperatures. [|I],0 It is generally expected that in the (Wigner) crystal 
state, electronic spins will be either ferromagnetically or antiferromagnetically ordered. [^,^ 
There has long been theoretical interest |^,^ in the rather subtle physics which determines 
how the electronic spins are ordered. For two-dimensional electrons in the Wigner crystal 
zero magnetic field a series iTIJHl of variational and Green's function Monte Carlo 
calculations have not led to definitive conclusions concerning the nature of the magnetic 
order. The energetically preferred spin ordering has been shown to depend very much on 
the lattice structure of the electron crystal and, unfortunately, the difference in energy 
between states with different spin order on a hexagonal lattice (which is expected to be 
the ground-state lattice of the Wigner crystal) is smaller than the accuracy of the Monte 
Carlo calculations §]. To our knowledge, there are no previous numerical studies of the 
spin structure of the Wigner crystal in the strong field regime. The variational Monte 
Carlo calculations cited above considered, in the strong-field limit, only the spin-polarized 
hexagonal lattice and investigated exchange, correlation and Landau-level-mixing effects [§. 

In this paper we discuss magnetic order for two-dimensional electrons in the limit of 
strong perpendicular magnetic fields where all electrons are confined to the lowest quantized 
kinetic energy Landau level. In this limit, the state of the electrons depends on the Landau 
level filling factor v rather than the electron density and, except for a narrow interval 
surrounding v = 0.2, the electrons form a Wigner crystal state for u smaller than ^ 0.23. 
[u = N/N^ where is the number of electrons and A^,^ = SB/^q = S/ (27r£^) is the Landau 
level degeneracy. Here S is the area of the system, B is the magnetic field strength, i is 
the magnetic length and $o = hc/e is the magnetic fiux quantum.) We find, partly on 
the basis of Hartree-Fock approximation (HFA) calculations, that in this regime the Wigner 
crystal state will always be ferromagnetic. The Hartree-Fock ground-state wave function 
does not contain the important correlations that give rise to the magneto-phonon and spin 
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waves modes of the crystal. We can check, however, the stabihty of the spin-polarized 
lattice by evaluating the magnon spectrum of the Wigner crystal using a time- dependent 
Hartree-Fock approximation (TDHFA). Within the limits of our numerical approach, we 
find that the polarized lattices remain stable at small filling factors. Moreover, in this limit, 
the spin-wave modes are very well described by a Heisenberg model where electrons are 
localized on their lattice site with an effective exchange integral Jtdhfa that we compute 
for different filling factors. From the value of this effective exchange integral, we can derive 
the spin-stiffness of the Wigner crystal state and compare its value with the spin-stiffness 
of ferromagnetic electron fluid states at nearby filling factors. This comparison shows that 
the spin-stiffness of the liquid is much larger than that of the solid. This is so because, the 
exchange energy is larger when the electrons are free to move around and come closer to 
each other. 

Our paper is organized as follows. In Section II we outline the formalism we use to 
perform the HFA and TDHFA calculations calculations for respectively the ground state and 
spin waves of the system. We are able to enormously simplify the calculations by adapting 
an approach we developed previously to the situation of interest here. In Section III we 
present numerical results for the magnetic ground state of the square and triangular Wigner 
crystal states and discuss differences between magnetic ordering tendencies in zero-field and 
strong-field limits. Our results for magnon dispersion relation are presented and discussed 
in Section IV. Section V contains a brief summary of this work. 

II. HARTREE-FOCK AND TIME-DEPENDENT HARTREE-FOCK 

APPROXIMATIONS 

A. Hartree-Fock Approximation Ground State 

The formalism outlined in this section is a straightforward generalization of one which 
we developed originally ||^ to describe phonon modes in the Wigner crystal state and which 
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has previously been generalized in other directions [|TO|,|TT| to describe double-layer quantum 
Hall systems and edge excitations of the Wigner crystal. We outline the main steps in the 
development of the formalism and refer the reader to Ref. for further details. 

We consider a two-dimensional electron gas (2DEG) in a magnetic field B = —Bz which 
is assumed to be strong enough so that we can make the usual approximation of considering 
only the lowest Landau level. In the Landau gauge, the Hamiltonian of the 2DEG is then 
(we set ^ = 1 throughout this paper) 

H = T.^-^xCa,x + ^Y: E E^(q)(^l|exp(^q■r)|X4) 

a,X ^'^ q Xi,...,X4 a,f3 

X (X2I exp(-iq ■ r) 1X3)0!^^^ 4^^2C/3,x3Ca,X4, (1) 
where = -|-(up), — (down) are spin indices and the lowest Landau level has the energy 

As usual, Uc = eB/m*c is the cyclotron frequency and m* and g* are the effective mass and 
(7-factor of the electron appropriate to the two-dimensional electron layer. For a finite 
system, the allowed values of the quantum number X are separated by 2ni'^/Ly. Neglecting 
the finite thickness of the two-dimensional electron layer, we take ^(q) = 27re^/g, the two- 
dimensional Fourier transform of the Coulomb potential. 

Making the usual Hartree-Fock pairing of the second-quantized operators in the Hamil- 
tonian of Eq. (|l]) and allowing for the possibility of broken translational symmetry and spin 
magnetic order in the ground-state, we obtain 

H = N^Y. e«Pa,a(o) + iv<^ E E Vo.A^)po.A^), (3) 

" q a,(i 

where we have introduced the operators 

Pa,/3(q) = E {-^xX - ci^xCf3,X+qyP, (4) 

X 

which are related to the density and spin operators by the relations 
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n{q) = iV^e-'?^^^/4[p++(q)+p__(q)], 
S\q) = iV<^e-^^^^/V+4q), 

The matrix elements of the Hartree-Fock self-consistent field in Eq.(^ are given by 

V^++(q) = [/7(q)-X(q)](p++(-q))+i/(q)(p__(-q)), 
Vl_(q) = [H{cO - X(q)] (p--(-q)) + H{q) (p++(-q)) , 

(6) 

V+.{q) = -X(q)(p_+(-q)), 
\/_+(q) = -X(q) (p+4-q)) , 
with the Hartree (if) and Fock (X) interactions defined by 

where /o(a;) is the modified Bessel function of the first kind and the factor (1 — (5q^o) comes 
from the neutralizing positive background. 

The ordered state is defined by the set of order parameters {(Pa,/3(q))} • In the case of 
interest here, i.e. for a Wigner lattice, these parameters are non-zero only when q = G, 
a reciprocal lattice vector of the crystal. To calculate the {pa^piq)) 's, we define the 2x2 
single-particle Green's function 

G^AX,X',t) = -(Tc„,x(r)4^,(0)), (8) 



and its Fourier transform Ga^^(q, r) by 

G„,^(q, r) = iV-i GaA^, X', r) exp [-^zg.(X + X')] Sx',x-,yP, (9) 
X,X' ^ -' 

so that 

{Pa,M)=Gf3A^,T = 0-). (10) 

Using the Heisenberg equation of motion ^(. . .) = [H—ixN, (...)] where p is the chemical 
potential of the electrons which we measure with respect to the kinetic energy of the first 



Landau level, we obtain the equation of motion for the single-particle Green's function (in 
an obvious matrix notation) 



{iUn + /i) / — A 



^ 0^ 



-tq X q 
2 ^ ^ 



V(q'-q)G(q',a;„)=/5q,o, 

(11) 



where Un is a fermionic Matsubara frequency, A = g*fiijB/2 and / is the 2x2 unit matrix. 

Eq. (^) is very general. For example, it can be used to consider complex spin-texture 
states such as the Skyrme crystal studied in Ref. [|10| where the average value of all three 
components of the average spin are space dependent. Although we concentrate, in this 
work, on simple spin-structure states where Eq. ([TT]) can be reduced to only one uncoupled 
equation, we explain here our numerical approach for the general case. 

We represent by qi, q2, qs, . . . qA? the wave vectors defining the ordered state (in principle, 
N oo but in the numerical calculation, a suitable cutoff is chosen for A^). We choose 
qi = and define the vector Ga,i3 = {Ga,p{qi) , GaA^2) , G^A^s) , ■ ■ ■ ,Ga,i3{^N))- Since in 

Eq. (|ll|), G++ (or G ) is coupled to G^+ (or G'+_ ) only, we can simplify Eq. ([lT|) by 

defining the 2A^— component vectors Gi = (G'++,G'_+) and G2 = iG+-,G ). We finally 

get a set of two coupled integral equations that we write in matrix form as 



-F 



(12) 



and 









\ 






-F 














- / 





(13) 



In these equations, I is the 2A^ x 2A^ unit matrix, 1 = (1, 0, 0, . . . , 0) and = (0, 0, 0, . . . , 0) 
are respectively the A^— component unit and nul vector, and F is the 2A^ x 2A^ matrix defined 
by 

(q,q') A+_(q,q') 



F = 



A(5q,q/ + 



^-+(q,q') 



-A5q,q, +A__(q,q') 



(14) 
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where 



^a,/3(q,q') = exp 



1 



-zq X q 



- q) 



(15) 



Note that since 74a,/3(q, q') = [A/3^Q,(q', q)]* , F is an hermitian matrix. It follows that 
Eqs. (|12|) and ( [T3| ) can be solved by making the unitary transformation F = UDU\ where 
UW = 1 and D is the diagonal matrix containing the eigenvalues of F. Following Ref. 
we have for the order parameters { i = 1,2, . . . , N ) 



(p++ 


(q.)) 




Ui,kUl^k, 




(q.)) 


V^^ — ^max 

— 2^k=l 


Ui+N,kUi^k-> 




(q.)) 


\ — ^max 

~ 2^k=l 




{p- 


(q.)) 


\ — ^max 

~ 2^k=l 


Ui+N,kUN+l^k 



(16) 



The value of /Cmax is obtained from the conditions 

(p++(0)) = u+ 

and 



(17) 



(18) 



the filling factors for spin up and down. It is easy to show from Eq. (|I6|) that, at T = QK, 
the following sum rules hold 



E[l(p++(q))l' + Kp+-(q))r 



(19) 



and 



E[l(p-+(q))l' + |(p--(q))r 



(20) 



We note that, except for simple cases such as the fully polarized or unpolarized crys- 
tals, the filling factors and z/_ are not known from the beginning. The only boundary 
conditions are the constraints 
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(p++(0)) + (p__(0)) = ^, 



(21) 



and 



(p+-(0)) = (P-+(0)) = 0, 



if A ^ 0. 



(22) 



Also, by definition, 



(p+-(q)) = (p-+(-q))*- 



(23) 



To find z/_|_ , z/_, Eqs. (|T2|,|TB|) must be solved self-consistently for a given value of z/+ and z/_ 
until a convergent solution is obtained. The process has to be repeated for different sets of 
and z/_ values until the lowest-energy solution is found. In this way we can determine 
the lowest-energy single Slater determinant consistent with any assumed translational and 
magnetic symmetry. The Hartree-Fock energy per particle of a particular ground-state 
configuration (with respect to the kinetic energy of the lowest Landau level) is: 

1 



E = A 



2u 



5:{[i7(q)-X(q)] [|(p++(q))P + l(P-(q))P] (24) 
q 

+i^(q) [(p++(q))(P-(-q)) + h.c] - 2X(q)|(p+_(q))r} 
In the case of a fully spin-polarized Wigner crystal, only {(p+_|_(G))} 7^ (G is a recip- 



rocal lattice vector) and Eq. (11) simplifies to 



[lu^ + fi-A) G'++(G, .;„) - F(G, G')G'++(G', a;„) = 5g,o, 



(25) 



where F is now the N x N matrix 



F(G, G') = exp 



1 



-iG X G 



'/J2 



V^^(G-G') 



with 



(26) 



V^^{G) = [H{G)-X{G)] (p++(G)) 



(27) 



The ground-state energy per particle, in this case, is simply 

= - A + ^ E [H{G) - X{G)] I (p++(G)) r 



(28) 
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B. Time-Dependent Hartree-Fock Approximation Collective Excitations 

To determine the collective excitation energies of the ordered state, we define the response 
functions 

Xap^si^l, q'; r) = -g{Tpaf3{(i, r)p^5(-q', 0)), (29) 

where pap = Paf3 — (Pafs)- By making use of the commutation relation of the operators 

Pa(3{^} and of the HF Hamiltonian of Eq. we obtain an equation of motion for the 

response functions that corresponds to the HFA which we denote by We get (repeated 
spin indices are summed over) 

[iQn + (Ca - e/3)] X%-ysi'i, ^n) = 

Sp,,e-^'^^^^''' (p„,(q - q')) - 5^,se''^^-^''' (p,^(q - q')) 

q" 

+ E V(q" - q)e^^^^^"''xL,.(q", q', ^0, (so) 

q" 

where fin is a boson frequency. 

To calculate the response functions in the Time-Dependent Hartree-Fock Approximation 
(TDHFA), and so include the correlations that give rise to phonons and magnons, we need 
to sum a set of ladder and bubble diagrams . The final equation for x can be expressed 
solely in terms the order parameters of the crystal phase! 

Xaisjsiq., q'; ^n) = Xa/sjsi'i, q'; ^n) 

+ E XafsU^, q"; fiO^(q")x5€75(q", q'; ^n), (3i) 

q" 

where the irreducible response function is given by 

(q> q'; ^n) = x%js (q, q'; ^n) (32) 
- E Xa/3Ke(q> q"; ^n)x{q")x^^^siq!', q ; fin)- 
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The spin and density response functions are obtained, as usual, from the analytic contin- 
uation iD,n cu + i5. The dispersion relation of the collective modes are then found by 
tracking the poles of the response functions at different values of the wavevector q in the 
Brillouin zone. 

In the case of a fully-polarized state, the only non-zero response functions are 

Xh \-,X-++- and X++++ and so the usual spin flip and density-density response functions 

are given by 



X+-(q, q'; 0„) = ^e-'' ^ /'e-^'^^^/^x+-+(q, q'; ^n), 



(33) 



and 



X^^(q,q';r)n) = ^X""(q,q';^^n) 



(34) 



They obey the TDHFA equations of motion 

E [^^n5q,q" - Ca (q, q") - Da (q, q") [i/(q") - X(q")]] X++++(q", q'; ^n) = /^^(q, q'), 

q" 

(35) 

and 

E mn - 2A)5q,q. - Cb (q, q") + Db (q, q") X(q")] X+-+(q", q'; ^n) = ^^^(q, q'), (36) 

q" 

where we have defined 



^A(q,q') = -2^sin [(q x q')£V2] , 
^B(q,q') = (p++(q-q'))e-^^'^''^'^''/', 



(37) 
(38) 



C^(q, q') = 2z (p++(q - q')) [H{ci - q') - X{q - q')] sin (q x q')£72 



(39) 



Cb{ci, q') = (p++(q - q')) X{q - q') cos [(q x q')^72] 

+ (p++(q - q')) [2iH{(i - q') - zX{q - q')] sin [(q x q')£V2 



(40) 
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(For a Wigner crystal, q — > k + G, q' k + G' etc. where k is a vector restricted to the first 
Brillouin zone of the crystal.) The problem of calculating the spin- flip and density-density 
response functions is then reduced to a matrix-diagonalization problem. The two response 
functions decouple. The matrix eigenvalues are the collective excitations associated with 
the two response functions, phonons in the case of x""" and magnons in the case of ■ 



III. HARTREE FOCK APPROXIMATION FOR THE GROUND STATE 

We first apply the above formalism to examine the nature of the magnetic order in 
the Wigner crystal ground state. In the Hartree-Fock approximation the ground state at 



strong magnetic fields always has broken translational symmetry. [|T^ This result of the 
Hartree-Fock approximation is an artifact. As we mentioned in the introduction, the true 
ground state has broken translational symmetry only for v < 0.23. Nevertheless, as 
we discuss further below the Hartree-Fock approximation does describes the ground state, 
reasonably accurately when the ground state is a Wigner crystal. Of course the Hartree- 
Fock approximation completely misrepresents the excitation spectrum of the Wigner crystal, 
since it misses the phonon and magnon collective modes captured by the time- dependent 
Hartree-Fock approximation. 



In two-dimensions. Coulomb interactions favor a triangular lattice for the Wigner 
crystal. We will find that the energy scale associated with magnetic order is much smaller 
than the Coulomb energy scale. We therefore expect the structure of the Wigner crystal to 
be the triangular lattice structure dictated by Coulomb interactions. We also expect that the 
interactions between the spins on the triangular lattice sites will be predominantly nearest 
neighbor since the overlap between wave functions on different sites is quite small in a strong 
magnetic field. We check this approximation below by comparing the dispersion relation of 
the spin waves in the TDHFA with that given by a Heisenberg model with only nearest- 
neighbour exchange coupling. The ground state for two-dimensional spin 1/2 particles with 
nearest-neighbour interactions on a triangular lattice is expected to have long range order for 
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both ferromagnetic and antiferromagnetic interactions. [|T5[ However, because of frustration, 
the order is rather subtle for the antiferromagnetic case. (The triangular lattice is not a 
bipartite lattice.) Our primary objective in this subsection is to determine whether the 
interactions is ferromagnetic or antiferromagnetic by comparing the energy of these two 
states. For that purpose it is more useful to consider the case of two-dimensional electrons 
on a square lattice since it is bipartite and both antiferromagnetic and ferromagnetic states 
have a simple structure. We do so even though the ground state of the two-dimensional 
electron solid does not occur in this structure. 



A. Maki-Zotos Wavefunction 



It is instructive to begin by generalizing the wavefunction for spinless electrons employed 



by Maki and Zotos |jl6| in their study of the strong field Wigner crystal. We define 

v]/ = (iV!)-i/Met K(r.)xL, 



(41) 



Here R, is the j-th lattice vector. 



V'R(r) 



1 



Rl 



: exp 



2i{xRy - yR^ 



(42) 



V2^ \ y 

is the lowest Landau level wavefunction [|17| for an electron whose quantized cyclotron orbit 
is centered on R, and Xn = (cos(6'/2), sin(^^/2) exp(i0)) is a spinor oriented in the n = 
(sin(6') cos(0), sin(6') sin((/)), cos(6')) direction. In this wavefunction the cyclotron orbits of 
electrons near different lattice sites are uncorrelated and the electron spin orientation at a 
given lattice site is arbitrary. In the range of u where the ground state is a Wigner crystal it 



is an excellent approximation [0,113] to ignore the lack of orthogonality between cyclotron 
orbits centered at different lattice sites. Making this approximation, it is easy to derive an 
expression for 



E 



(^lEi<,e2| 


ri - 


-1 







(43) 
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In Eq. we have dropped the Zeeman energy which can easily be added if the electronic 
(7-factor is non-zero. The fact that the kinetic energy, taken as the zero of energy above, 
is quantized is important in determining the favored magnetic order. Following Maki and 
Zotos we find that 



^ = ^E /(|R,-R.|)-(^^^^)^Mz(|R,-R.|) 



1 



(44) 



where 



HR) = ( y ) ^ exp {-R^Sf) Jo {R'/Sf) , (45) 



and 



JMziR) = exp(-i?V4£2)/(i?), (46) 

In this equation I{R) and J{R) are respectively the direct and exchange two-body matrix 
elements of the Coulomb interaction for lowest-Landau-level cyclotron orbits whose centers 
are separated by R. The explicit expression for the matrix element in the exchange term, 
which is sensitive to the relative orientations of the spins on the two sites, is 

7 no 2f. /• , ,^R,WV'R,(r')^R,(r')^R.(r) 
Jmz (|Rj - = e' J dr J dr' — yyyj • ('^^) 

For i? > 

JMziR)-^exp{-RyAf). (48) 

This approximate expression for Jmz{R) is accurate to better than 5% even for neighboring 
sites over the range of Landau level filling factors where the ground state is a Wigner crystal. 

It is evident from Eq.(pD that if the ground state is approximated by the Maki-Zotos 
wave function, a ferromagnetic state in which all spins are parallel will be energetically fa- 
vored. The energy increase when the relative orientation of spins on two sites separated by 
R changes from parallel to antiparallel is J{R). For similar single-Slater-determinant varia- 
tional wave functions at zero magnetic field, the tendency would be to favor antiferromag- 
netic orientations on neighboring sites [^-§1 except possibly when multi-site ring exchanges 
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become important. (Multi-site ring exchanges are less important for the strong magnetic 
field Wigner crystal because magnetic confinement results in orbitals which are more strongly 
localized around lattice sites.) In the weak field case, having opposite spins on neighboring 
sites reduces the kinetic energy density required by the Pauli exclusion principle in the region 
between the sites. In the strong magnetic field limit, the kinetic energy is quantized and 
is independent of the spin-configuration so this mechanism favoring antiferromagnetism is 
not operative. Nevertheless the Maki-Zotos wavefunction is a single-Slater-determinant and 
conclusions based upon its use should be examined critically. It is known, for example, that 
correlations can result in spin-singlet fiuid ground states |jl9| whereas the Hartree-Fock ap- 



proximation would always predict ferromagnetic ground states. At zero magnetic field, the 
contribution from low-energy, long-wavelength phonon modes to the zero-point motion gives 
rise to long-range correlations which, for example, make the static structure factor vanish 
more quickly (oc g^/^) than it would for a system with short range interactions. At strong 
magnetic field, even stronger correlations which make the static structure factor vanish as 
result from the contribution to the zero-point motion of the collective cyclotron mode of 
all electrons. (In a Jastrow-Slater variational wavefunction such as that used by Zhu and 
Louie [^] the correlation factors would have to have a logarithmic spatial dependence in 
order to capture the correct long-distance ground-state correlations.) We cannot completely 
rule out on the basis of our calculations the possibility that correlations could invalidate our 
conclusion that the ground state is ferromagnetic. However we consider this to be extremely 
unlikely. 

B. Self- Consistent Hartree-Fock Calculations 

One possible mechanism in favor of antiferromagnetism is the possibility of spreading 
the charge associated with a given lattice site more widely in the case of antiferromagnetic 
configurations which could reduce the electrostatic energy. To probe the competition a little 
more deeply we have performed self-consistent Hartree-Fock calculations, based on the for- 
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malism of the previous section, comparing the energy of ferromagnetic and antiferromagnetic 
states on a square lattice. We now discuss the results of these calculations. 

In the Hartree-Fock approximation, the spin-order is unidirectional on a square lattice 
for both ferromagnetic and antiferromagnetic interactions. We choose a spin-quantization 
axis which is along the direction of the Zeeman coupling if one is present and is otherwise 
arbitrary. This allows us to set the order parameters which are off-diagonal in the spin indices 
to zero and simplify our calculation. Let oq be the lattice constant of the ferromagnetic 
square lattice with density n = I/oq such that 2'imC.'^ = v. In the antiferromagnetic case 
we assume that the spin density is oppositely directed on the two sublattices (which have 
lattice constant v^oq, and have a relative shift of a =^/2 ^ Oq) so that 

(p__(G)) = e-'^-(p+4G)) (49) 

where G is a sublattice reciprocal lattice vector (with modulus |G| = 27r/ v^ao)- (We choose 
our coordinate system so that the primitive lattice vectors of the sublattice are along the 
Cartesian axes.) Eq. (|TT|) can again be simplified to a single equation 

(2cu„ + p) G++(G, un) - F(G, G')G++(G', Un) = 5g,o, (50) 

G' 

where 

1 



F(G, G') = exp 



-iG X G'e^ 

2 



(l + e-«-) H{G) - X(G)] (p++(-G)) • (51) 
The ground-state energy per particle becomes 

E^- = W WG) (1 + cos(G - a)) - X(G)] | (p++(G)) ^ • (52) 

^ G 

We have solved these equations self-consistently. Because of the variational nature of the 
Hartree-Fock approximation, these solutions provide us with the lowest energy single-Slater 
determinant consistent with the assummed magnetic and translational broken symmetry. 
In particular, the solutions to these equations will always give a lower energy than the 
energy for the corresponding Maki-Zotos wavefunction. The optimization process implicit 
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in obtaining a self-consistent solution of the Hartree-Fock equations results in cyclotron 
orbits on each lattice site which are distorted by their average environments, including their 
magnetic environments, in a way which minimizes the total interaction energy. It is still 
true, however, that the cyclotron orbits on different sites are not correlated with each other. 
The error introduced as a consequence can be estimated by using a harmonic approximation 
for the strong field Wigner crystal, which is reasonably accurate from an energetic point of 
view throughout the regime where the ground state is an electron crystal. In the harmonic 
approximation the many-body Schrodinger equation can be solved exactly and the ground- 
state energy is the sum of the classical Madelung energy and the quantum zero-point energy 
i.e. 

Eharmonic = -0.78213z/^/2 ^ 0.2410 1//=^/^^ (53) 

for the hexagonal lattice. The Hartree-Fock approximation describes the Madelung term 
exactly (in the limit u 0, the HFA energy coincides with the classical energy of a point 
lattice) and overestimates |^| the zero-point energy by approximately 25% (at u = 0.2 ). 



The results of our calculations are summarized in Tables land 2. In Table 1, we list 
the ground-state energy per electron in the HFA for the square lattice antiferromagnetic 
(SLA) and ferromagnetic (SLF) states as well as for the triangular lattice ferromagnetic 
state (TLF). Table 2 contains a similar calculation using a simplified form of the Maki and 
Zotos wave function where we have neglected the overlapping between two wave functions 
centered on different sites so that the single-electron density can be approximated by 

Hr)) = E \Mr)f = ^,T.e-''-''^''^''\ (54) 

i i 

or, equivalently, in the ferromagnetic case 

{piG)),,, = V e-^'''/\ (55) 

(For the SLA case, v v/2 and the G's are replaced by the sublattices reciprocal lat- 
tice vectors.) We use the order parameters defined by Eq.(^) in Eg . (pSf) to compute the 

16 



Maki-Zotos ground-state energies tabulated in Table II. We remark that this procedure is 
exactly equivalent to computing Eq. ( PI) (when the interaction with a positive homogeneous 
background of charges is added to this last equation). Note also that these results include 
only the Coulomb energy. These Tables report also the results of calculations performed for 
filling factors where the ground state is not believed to be a Wigner crystal. These large 
V results are intended to illustrate trends in the self-consistent Hartree-Fock approximation 



solutions and not to be physically realistic. ^1 



We see that for the larger filling factors, the difference in energy between Hartree-Fock 
square lattice ferromagnetic and antiferromagnetic states, AEgpin = {Eslf — Esla) / Eslf 
agrees quite closely with what would be predicted by the Maki-Zotos wavefunction (Eq. 
(^)). At smaller filling factors, however, the Hartree-Fock energy difference between these 
two spin states is much bigger than what would be predicted by the Maki-Zotos wave- 
function. The energy reduction due to the added variational freedom compared to the 
Maki-Zotos wavefunction is larger for the ferromagnetic state than for the antiferromagnetic 
state and this leads to an increased energy difference between the two states. (See also 
Table 3 where, as discussed below, Jmz is proportional to {Eslf — Esla) evaluated with 
the Maki-Zotos wavefunction and Jhfa is proportional to {Eslf — Esla) evaluated in the 



HFA.) ||22| In both the Makis-Zotos and HF approximations, the ferromagnetic state has the 
lowest energy. In Fig. 1, we plot the (i2j^erence in density: (^(r))^:^^^— (r7,(r))^y^^ for the SLF 
and SLA at filling factor u = 1/8. It is clear from this figure that the HFA minimizes the 
Coulomb energy in both the SLF and SLA cases by removing charges along the direction of 
the nearest-neighbor sites and putting them along the direction of the next-nearest-neighbor 
sites. This is just what we expect at such a small filling factor where overlap between wave 
functions on different sites is very small and the ground-state energy is dominated by the 
(direct) Coulomb interaction. According to our calculations more charge redistribution oc- 
curs in the ferromagnetic case. For the triangular lattice, a similar calculation gives a much 
smaller difference in densities reflecting a loss in the variational freedom due to the higher 
coordination number of the triangular lattice. 
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We remark that similar self-consistent calculations for a single-band Hubbard model at 
half-filling would find the antiferromagnetic state to be lower in energy, correctly reflecting 



the superexchange coupling in that system. [^] We also see that, as anticipated above, the 
difference between the square lattice ferromagnetic state and the triangular lattice ferro- 
magnetic state energy AEcouiomb = {Etlf — Eslf)/ Etlf is much larger than the difference 
between ferromagnetic and antiferromagnetic states on the same lattice. This energy differ- 
ence is almost constant over the range of filling factors considered here, decreasing slowly as 



u decreases. For z/ ^ 0, AEcouiomb approaches its Madelung energy value, 12 0. 



0. 



IV. COLLECTIVE MODE CALCULATIONS 

As we described above, the ground-state order parameters can be used to calculate the 
spin-wave collective modes. In the ferromagnetic ground state, we showed that the density 
response function, x""', or equivalently the longitudinal spin response function, = X""/4, 
are uncoupled from the transverse spin response function . The poles of the longitudinal 
spin response function are nothing but the phonons of the Wigner crystal for which we have 
already computed the dispersion relation in Ref. The transverse spin excitations are 
the magnon collective modes of the ferromagnetic Wigner crystal. Figs. 1 and 2 show the 
TDHFA dispersion relations for the SF and TF lattices, respectively, at different values of 
the filling factor. (When the electronic g-factor is non-zero all magnon collective modes 
energies are increased by the Zeeman gap g*fiBB. In the simple colinear states (SLF, SLA 
and TLF) that we consider here, a small Zeeman term has no effect on the calculation of 
the order parameters.) We remark that, because of the numerical approach used in this 
work, it is not possible to obtain the dispersion relations at very small filling factors without 
having to consider a prohibitively large number of reciprocal lattice vectors. For the square 
lattice, we were not able to obtain accurate results for u < 1/5 while for the triangular 
lattice u = 1/7 was the lower limit. 

It is interesting to compare these magnon dispersion relations with the dispersion relation 
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of the spin waves of the Heisenberg model where the spins are locahzed on the lattice sites 
and the Hamiltonian is given by 

H = -JJ2S,-S,+s. (56) 

iS 

(The summation is over the lattices sites i and the nearest neighbours of the lattice. Note 
that this convention for the exchange constant results in double-counting each neighbour 
pair.) In that case, 

a;(k) =2Jz/os(l-7k), (57) 

with 

7k = - E (58) 

<5 

and s = 1/2. 

The solid lines in Figs. 1 and 2, show the dispersion relation obtained from a nearest- 
neighbor interaction Heisenberg model with the interaction strength chosen to reproduce the 
TDHFA numerical results. The fit is quite good and becomes almost perfect at smaller filling 
factors. (The discrepancy at z/ = 1/3 can be improved by fitting with non-zero next-nearest- 
neighbour coupling.) The exchange integral Jtdhfa obtained, in this way, from the TDHFA 
dispersion relation, for different filling factors, is listed in Table 3 (square lattice) and Table 
4 (triangular lattice) . These tables also show values of the exchange integral calculated 
in two other ways. From Table 1, we see that the ground-state energy is minimal for the 
polarized square lattice. We can estimate the strength of the exchange coupling from this 
energy difference as follows. We assume that, as far as the magnetic degrees of freedom are 
concerned, the Hartree-Fock solution yields an Ising approximation to the antiferromagnetic 
ground state, i.e. it does note capture the quantum fluctuations which would be present in 



a true antiferromagnetic ground state. It is then easy to see [24] that, for the square lattice 



JuPA = ^^'^-^'^r (59) 
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This expression assumes that non-nearest neighbour exchange couphng is neghgible. (Note 
that a similar calculation is not possible for the triangular lattice because of frustration.) 
We can also compute the exchange integral directly from the Maki-Zotos wavefunction ex- 
pression, Eq.(^). This gives Eq.(^5D or, using 27m£^ = v with n = l/aa^ (a = 1 (SF) or 



Note that if the Maki-Zotos energies are used in Eq. ( ^9|) instead of the HFA energies, Jmz 
is recovered exactly as long as non-nearest-neighbour interactions are negligible. 

For the TLF, the exchange integral obtained from the TDHFA is only slightly smaller 
than Jmz- For the square lattice, the exchange integral obtained from the TDHFA is smaller 
than both Jmz and the HFA value over the range of u where we are able to complete 
calculations. At smaller filling factors, since the difference between the two spin states 
decreases faster with the Maki-Zotos wavefunction than with the HFA, the HF value of J is 
much larger than the Maki-Zotos value. We expect the TDHFA result to remain close to the 
HFA result in this regime. Unfortunately, we cannot check this assumption numerically since 
we cannot compute the TDHFA value of J at smaller filling factor (the matrix size becomes 
prohibitively large). However, since the TDHFA is obtained from a functional differentiation 
of the HFA, our assumption seems reasonable. In any case, the present result shows clearly 
that, in the Wigner crystal, the interaction between spins, at small filling factor, is mainly 
from nearest-neighbors. 

In the small wave vector limit, the Heisenberg dispersion relation on the triangular 
lattice is given by uj{k) = lJ{kaof = D{Hf where D = l-n^flJlv in units of e^ji. With 
J given by Jtdhfa as calculated above, we find that D = 0.024 at z/ = 1/3 and D = 9.3 
X 10~^ at u = 1/5. These values should be compared with those of the liquid state where 
D = Ani'^ps/i' where ps is the spin-stiffness. For the liquid state it is possible to express the 
spin-stiffness in terms of the pair correlation function |TP| , |25[| and this has been evaluated 
using an hypernetted-chain-approximation for the liquid state pair-correlation function in 



a = 



73/2 (TF)) 




(60) 
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Ref. For the liquid state we find that, D = 0.035(eV^) at z/ = 1/3 and D = 0.015(eV^) 
at u = 1/5. We see that the spin-stiffness is larger for the liquid state and increasingly so as 
the filling factor decreases. This result is consistent with the view of the strongly correlated 
electron states as quantum melted crystals of electrons whose size is smeared on a magnetic 
length scale by rapid cyclotron motion. When long-range order is lost, the cyclotron orbits 
will overlap more strongly on average and the relative spin-orientation of nearby electrons 
will assume a larger importance. 

V. SUMMARY 

In the strong- magnetic-field limit, we have argued that the lowest-energy spin state of 
the Wigner crystal is the ferromagnetic state. Our conclusion is based in part on a compar- 
ison of ferromagnetic and antiferromagnetic state energies of square-lattice Wigner crystal 
states calculated in the Hartree-Fock approximation and in part on the observation that the 
superexchange mechanism, which tends to favor antiferromagnetism, is absent in the strong 
magnetic field limit. The spin-wave dispersion relations, which we compute in the TDHFA, 
show that the ferromagnetic lattice is stable at filling factors where crystallization occurs. In 
this limit, the interactions between spins on the lattice are dominated by nearest-neighbor 
exchange coupling. Our results appear to show that small distortions of the wavefunctions 
for electrons on one lattice site, due to their interactions with electrons on nearby lattice 
sites, are responsible for a large relative increase in the small exchange couplings at small 
filling factors. A comparison with the liquid state at filling factors z/ = 1/3 and u = 1/5 
shows that the spin-stiffness of the ferromagnetic liquid states which occur at these filling is 
substantially larger than that of corresponding crystal states. In closing we remark that the 
recent sucessful application of nuclear-magnetic-resonance methods to two-dimensional 
electron system, suggests that the magnetic properties of two-dimensional electron systems 
in the regime where the Wigner crystal state occurs will soon be open to experimental in- 
vestigation so that our conclusions can be tested. These experiments should open up a host 



21 



of interesting new questions, related to disorder and thermal fluctuations. 
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TABLE CAPTIONS 



Table 1 Energy of the square and triangular ferromagnetic states (SLF,TLF) and of the 
antiferromagnetic state on the square lattice (SLA) in the Hartree-Fock approximation 
in units of /^. The relative energy difference between the ferromagnetic and antifer- 
romagnetic states on the square lattice is given by AEgpin while the relative difference 
in energy between the ferromagnetic state on the triangular and square lattice is given 

by A.E Coulomb ■ 

Table 2 Energy of the square and triangular ferromagnetic states (SLF,TLF) and of the 
antiferromagnetic state on the square lattice (SLA) in the Maki-Zotos approximation 
in units of e^/£. The relative energy difference between the ferromagnetic and antifer- 
romagnetic states on the square lattice is given by AEspm while the relative difference 
in energy between the ferromagnetic state on the triangular and square lattice is given 

by AEcoulomb- 

Table 3 Value of the exchange integral, in units of e^/£, and on the square lattice obtained 
from various approximations: Jmz is from the definition of the exchange integral given 
in Eq. ([60|), Jhfa is obtained from the energy difference between the ferromagnetic 
and antiferromagnetic states (Eq. (pOj)), and Jtdhfa is obtained by fitting the TDHFA 
spin-wave dispersion relation with the spin-wave dispersion relation of the Heisenberg 
model (see text). 

Table 4 Value of the exchange integral , in units of e^/£, and on the triangular lattice 
obtained from the various approximations described in Table 3. 
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FIGURE CAPTIONS 



Fig. 1 Difference in densities between tfie Hartree-Fock and Maki-Zotos approximations 
at filling factor v — 1/% and in units of I/oq for (a) the ferromagnetic square lattice 
and (b) antifcrromagnetic square lattice. For the ferromagnetic state, the lattice sites 
are indicated by gray circles. For the antifcrromagnetic state, sites on one sublattice 
are indicated by black circles and sites on the other sublattice are indicated by empty 
squares. Note that the orientation of the lattice differs in (a) and (b). 

Fig. 2 Dispersion relation of the spin waves of the ferromagnetic square lattice obtained from 
the TDHFA and from the Heisenberg model at different filling factors. The dispersion 
relation is plotted along the edges of the irreducible Brillouin zone of the square lattice 
with lattice spacing a. In units of 27r/a, F = (0, 0), J = (1/2, 1/2), X = (1/2, 0). 

Fig. 3 Dispersion relation of the spin waves of the ferromagnetic triangular lattice obtained 
from the TDHFA and from the Heisenberg model at different filling factors. The 
dispersion relation is plotted along the edges of the irreducible Brillouin zone of the 
triangular lattice with lattice spacing a. In units of 27r/a, F = (0, 0), J = (1/3, 1/^3), 
X = (l/V3,0). 
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p 




EsLA 




^E^i„ii 


-^Ec,„ii„i„i, (%) 


1/3 


-.3857672 


-.3823134 


-.3884928 


.90 


.70 


1/4 


-.3484399 


-.3478963 


-.3511452 


.16 


.77 


1/5 


-.3196321 


-.3195124 


-.3219969 


.037 


.73 


1/6 


-.2964916 


-.2964471 


-.2985717 


.015 


.70 


1/7 


-.2775104 


-.2774874 


-.2793787 


.0083 


.67 


1/8 


-.2616473 


-.2616346 


-.2633555 


.0049 


.65 


1/10 


-.2365415 


-.2365379 


-.2380218 


.0015 


.62 



R. Cote and A.H. MacDonald, Table 1 
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p 


EsLi- 


EsLA 




^E^i„ii 


-^Ec,„ii„i„i, (%) 


1/3 


-.3858614 


-.3814086 


-.3885208 


1.15 


.68 


1/4 


-.3482640 


-.3474804 


-.3511413 


.23 


.82 


1/5 


-.3194893 


-.3193454 


-.3219900 


.045 


.78 


1/6 


-.2964083 


-.2963812 


-.2985680 


.0091 


.72 


1/7 


-.2774646 


-.2774594 


-.2793770 


.0019 


.68 


1/8 


-.2616223 


-.2616213 


-.2633548 


.00038 


.66 


1/10 


-.2365341 


-.2365340 


-.2380216 


.000016 


.63 



R. Cote and A.H. MacDonald, Table 2 



29 



;/ 


■/uz 


■hi FA 


JroiiFA 


1/3 


.22 X 10-2 


.17 X 10-2 


.13 X 10-2 


1/4 


.39 X 10-3 


.27 X 10-3 


.24 X 10-3 


1/5 


.72 X 10-^ 


.60 X 10-^ 


.43 X 10-^ 


1/6 


.14 X 10-^ 


.22 X 10-^ 




1/7 


.26 X 10-5 


.11 X 10-^ 




1/8 


.50 X 10-6 


.63 X 10-5 




1/10 


.19 X 10-^ 


.18 X 10-5 





R. Cote and A.H. MacDonald, Table 3 
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■/uz 


Jtdiifa 


1/3 


.99 X 10-^ 


.74 X 10-=^ 


1/4 


.14 X 10-3 


.12 X 10-3 


1/5 


.20 X 10-^ 


.17 X 10-^ 


1/6 


.29 X 10-5 


.25 X 10-5 


1/7 


.44 X 10-6 


.40 X 10-6 


1/8 


.67 X 10-^ 




1/10 


.16 X 10-^ 





R. Cote and A.H. MacDonald, Table 4 
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R. Cote and A.H. Macdonald, Fig. 1(a) 
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